Effect of intracellular diffusion on current— voltage curves in potassium channels 
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We study the effect of intracellular ion diffusion on ionic currents permeating through the cell 
membrane. Ion flux across the cell membrane is mediated by special proteins forming specific chan- 
nels. The structure of potassium channels have been widely studied in recent years with remarkable 
results: very precise measurements of the true current across a single channel are now available. 
Nevertheless, a complete understanding of the behavior of the channel is still lacking, though molec- 
ular dynamics and kinetic models have provided partial insights. In this paper we demonstrate, by 
analyzing the KcsA current-voltage currents via a suitable lattice model, that intracellular diffusion 
plays a crucial role in the permeation phenomenon. The interplay between the selectivity filter be- 
havior and the ion diffusion in the intracellular side allows a full explanation of the current-voltage 
curves. 
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I. INTRODUCTION 

Potassium currents across nerve membranes have been 
widely studied (see, e.g., [TH3] and the reviews |1HH])- 
Ionic channels selecting potassium currents are present 
in almost all types of cells in all organisms and they play 
many important and different functional roles. 

Different types of measurements ITUI provide a very 
detailed description of the behavior of potassium chan- 
nels. In general, less is known on their structure [6]. 
KcsA, a potassium channel from Streptomyces lividans, 
is the first ion channel whose structure has been identified 
via X-ray crystallography [llj . 

All ionic channels form selective pores in the cell mem- 
brane which open and close and, when in the open state, 
allow permeation of a selected ionic species (potassium 
in K+-channels). Their ability to open and close, i.e., 
gating, and their ability to allow the flux of a particular 
ionic species, i.e., selectivity, are not yet completely un- 
derstood. A lot is known in the case of the KcsA channel 
and, with some care, can be extrapolated to the whole 
family of K"*" channels. 

In KcsA, see for instance the detailed description in 
[SI I12j . gating is realized via four crossing transmem- 
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brane helices on the intracellular side. When in the open 
state, a spherical water-flUed cavity of diameter about 
10 A widens on the intracellular side of the channel up 
to the membrane plane. There a 10-15 A long and 
3 A large channel containing the selectivity filter con- 
nects the cavity to the extracellular side. 

When the channel flips to the open state a solution 
with the cytoplasm concentration reaches the entrance 
of the channel. Part of the ions permeates through the 
channel leaving an ion depleted region close to the pore. 
The typical time needed to restore the intracellular con- 
centrations close to the pore will depend strongly on the 
diffusion process of ions inside the cell. We can then 
imagine that the current flowing through the channel will 
depend both on the diffusion of ions in the cytoplasm and 
on the behavior of the selectivity filter. 

The problem of computing the permeation current in 
the open state, namely, the so called true current, has 
been approached theoretically by a large variety of meth- 
ods. Molecular dynamics studies [131 HI! gi'^^ ^ ^^^Y 
tailed description, but they usually do not provide macro- 
scopic currents estimates due to the too small involved 
time scale. Kinetic models [T5lf20] give very useful infor- 
mation, since electro-physiological experiments are per- 
formed over time scales much longer than the atomic one, 
with the drawback of the extreme simplification on the 
structure of the channel. 

In the recent literature it has been examined the pos- 
sibility to validate these kinetic models by comparing the 
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predicted behaviors with those observed experimentally. 
In particular the models have been tested against their 
ability to predict the current-voltage behavior. In this re- 
spect very accurate experimental results have been pub- 
lished for different potassium channels, see [H], [55], and 
[55] for the KcsA, MaxiK, and the Shaker channels, re- 
spectively. 

Models such as those quoted above describe to some 
extent the dynamics of ion permeation through the selec- 
tivity filter, the concentration of the ion in the cell is in- 
troduced in the model as a constant parameter. In other 
words in those studies the diffusion of the ions inside the 
cell is not taken into account. Our opinion is that dif- 
fusion, as explained above, must take an important part 
in the permeation phenomenon. In |24| , inspired by |20j , 
we introduced a model where the channel is lumped to 
a two state stochastic point system and the interaction 
between the dynamics of the ions inside the cell (diffu- 
sion) and that of the selectivity filter itself is taken into 
account. The channel is then thought of as part of the 
cell more than as an isolated structure. In that paper 
both an analytical and Monte Carlo study showed the 
possibility to achieve gating via selection. 

In this paper we examine the possibility to predict the 
behavior of the current -voltage curves (graph of the per- 
meation current vs. the external potential difference ap- 
plied to the membrane) on the basis of a similar model. 
A modification is needed to take into account the effect of 
an external voltage difference through the cell membrane. 
The model is thus defined to mimic the three effects that 
seem to be the most relevant in the process: (i) diffusion 
of the ions inside the cell; (ii) dynamics of the selectivity 
filter; (iii) dynamics of the ions inside the selectivity fil- 
ter. We compare the current -voltage behavior predicted 
by our model with experimental results from j21H23| and 
find a very good agreement both for the dependence of 
the current on the ion concentration and on the external 
voltage. 

Ion diffusion (item (i)) inside the cell is modeled as a 
symmetric random walk on a finite line. We use a one di- 
mensional system in order to compute all the interesting 
quantities explicitly [53] Appendix B]. As we shall com- 
ment later the dimensionality of the system affects only 
our estimate of the diffusion coefficients of ions inside the 
cytoplasm. Corrections will be introduced to compute 
the three-dimensional value of the diffusion coefficient. 

We remark that point (i) is the real distinguishing fea- 
ture of our model. As we shall discuss later, the introduc- 



tion of ion diffusion in the intracellular region and hence 
of the depletion phenomenon, will allow a full description 
of the permeation current behavior with respect to both 
external voltage and ion intracellular concentration. 

One of the two boundary points of the finite line where 
ions diffuse is reflecting, whereas the other mimics the 
selectivity fllter. The dynamics of the pore (item (ii)) 
is assumed to be stochastic. The pore jumps randomly 
between two states, the low and the high-affinity one. 
The dynamics of the pore is independent from that of 
the ions inside the cell. The chances that an ion has to 
enter the pore depend both on the ion species and on the 
pore state; in this way selection is implemented in the 
model. 

This description is quite faithful to the real behavior 
of the selectivity filter. Indeed, two possible states are 
possible for the filter [25j, the low and the high-affinity 
one. When the filter is in the low-affinity state perme- 
ation is favored, but no ionic species is preferred. Thus, 
in order to realize selectivity in an efficient way, the filter 
has to jump between these two states. To our knowledge 
the time fraction spent by the filter in the low affinity- 
state is not known experimentally; the value predicted 
by analyzing the current-voltage curves via our model 
depends on the ionic concentration in cytosol and is of 
order 10~^. This result seems to be coherent with the 
qualitative description given in |201 Figure 5]. 

The dynamics of the ions inside the filter (item (iii)) 
is not modeled in detail, we just assume that a particle 
inside the pore can either exit the system or reenter the 
cell with a fixed probability. This ejection probability is 
chosen as a function of the voltage difference across the 
cell membrane. Our modeling of the dynamics inside the 
filter is, thus, reduced to the choice of this function. It is 
worth noting that the model is able to reproduce accu- 
rately the experimental results with different choices of 
this function, that is to say with different descriptions of 
the dynamics of the ions inside the filter. More precisely 
we shall see that it is possible to explain the experimental 
results via different models for the dynamics of the ions 
inside the filter; when a different functional behavior for 
the ejection probability in terms of the external voltage is 
chosen, a different value of the typical time spent by the 
filter in the low-affinity state is found. In other words 
different models of the dynamics of the ions inside the fil- 
ter are allowed provided the low-affinity state probability 
is changed suitably. 

Indeed, we shall find very good results either by choos- 
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ing the ejection probability function according to very 
well known and studied theories )17j or by assuming that 
the ejection probability for a potassium ion trapped in- 
side the filter is a power law function of the applied (suit- 
ably renormalized) voltage difference. In both cases our 
model will predict curves fitting cleanly the experimental 
result for the "true" ion current with reasonable values of 
the intracellular diffusion coefficient; but different values 
of the low-afhnity state probability will be predicted. A 
more precise knowledge of this probability would enable 
us to discriminate among different mechanisms. 

Moreover, we remark that our model will be also able 
to explain the behavior of ionic species different from 
potassium by choosing properly the physical parameters 
appearing in the ejection probability function of the ex- 
ternal voltage difference. 

We finally note that some care has to be used when 
the theoretical results are compared with experiments. 
In the biological literature, see for instance [9], two dif- 
ferent types of current are reported, the apparent and 
the true one. This is connected with the peculiar behav- 
ior of ionic channels: two different states are observed, 
closed and open channel. In the first state a zero cur- 
rent is observed, in the second one the outgoing current 
fluctuates randomly on a very short time-scale and a not 
zero average current, called true current, is measured. 
The true current can be accessed experimentally only if 
the time resolution of the instruments is good enough 
to distinguish neatly between the open and the closed 
time intervals. When this is not the case a different cur- 
rent, called apparent current, is measured. This current 
is smaller than the true one since the instruments average 
the current also on time intervals in which the channel is 
closed. 

The model we propose here is intended to mimic the 
behavior of the ionic channel in the open state. The 
average outgoing flux that will be computed has then to 
be thought as a prediction for the true current. 

In Section |ll] we introduce the lattice model. Sec- 



tion III is devoted to the analysis of the experimental 
data. Our results will be discussed in Section [TV] where, 
in particular, we shall comment on the physical mean- 
ing of the fitting parameters introduced previously. Our 
conclusions are briefly summarized in Section [V] 



II. THE MODEL 

The intracellular region is modeled via a finite one- 
dimensional lattice A C Z with L sites. Two sites of 
Z are said to be nearest-neighbors if and only if their 
mutual Euclidean distance is equal to one. The bound- 
ary dA of A is the collection of the two sites of Z not 
belonging to A and neighboring one of the site of A. One 
of the two points of the boundary of A is called pore, and 
denoted by P. 

One ionic species performs an independent symmetric 
random walk on the lattice with reflectivity conditions on 
the site of i9A different from the pore. The particles on 
the site neighboring the pore behave in a special way that 
will be described below; due to the peculiarity of such a 
rule the walkers will turn out to be not independent. The 
number of walkers is denoted by N. 

The fact that the walkers are independent on A means 
that the position of a particle does not affect the motion 
of the others, in particular no constraint to the number of 
particles on each site is prescribed. The fact that the ran- 
dom walks are symmetric means that each jump between 
two neighboring sites of A is performed with probability 
1/2. Since we assumed that the boundary is reflecting, 
particles in the site neighboring dA \ {P} can stay in the 
same site with probability 1/2. 

Two states are allowed for the pore: high-affinity and 
low-affinity. The pore switches between the two states 
randomly; the probability that the pore is in the low- 
afhnity state is denoted by p e [0, 1]. Moreover the pore 
can be either free or occupied by an ion. The behavior of 
the particles on the site Q neighboring the pore depends 
on the state of the pore itself as it is precisely stated 
in figure [T] The idea is the following. If the pore is 
occupied by a particle, no other ion can enter it. If the 
pore is (free) in the low-affinity state, particles can enter 
it; when they enter the pore they immediately dissociate 
so that, with probability v G [0, 1] they exit the system 
while they reenter A with probability 1 — v. If the pore is 
free and in the high-affinity state, one particle can enter 
it, but, once entered, it remains there until the pore state 
changes to the low-affinity one. When this happens the 
ion dissociates with the same rule described above. As 
noted above, due to the pore rule, the walkers are not 
independent. 

Whenever an ion exits the system another particle is 
put at random with uniform probability 1 /L on one of the 
L sites in A so that the number of ions is kept constant. 
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The model described above is implemented with the 
Markov Chain described in detail in the appendix |X] An 
iteration or sweep of the chain is the collection of the 
steps that are performed at each time n. 

As we have explained in the Introduction, the aim of 
this paper is that of computing the outgoing ionic cur- 
rent. This quantity is related to the number of particles 
that exit the system. We let M{n) be the number of par- 
ticles which exited the system in the time interval [0,7i]. 
Moreover, we let the flux at time n be M{n)/n. Since 
we defined the stochastic process in such a way that the 
number of ions keeps constant in the volume A, M(n) /n 
approaches a constant value / for n ^ oo. This quantity, 
that we shall call the outgoing flux, is expected to be pro- 
portional to the "real" current measured experimentally. 

The existence of this limiting flux can be deduced by 
remarking that the chain is irreducible and that the space 
state is finite. So we have that there exists a stationary 
measure for the process and that, in the limit rt — )> oo, 
the time dependent quantities tend to the corresponding 
quantities averaged against the stationary measure. 

By exploiting one-dimensionality, the model can be 
solved analytically and the outgoing flux can be com- 
puted explicitly. In our model the pore is modeled in a 
very simple fashion, indeed it is just a two state Bernoulli 
process; the main difficulty is, obviously, the interaction 
between the random walk inside the volume A and the 
pore itself. We consider the stationary state reached by 
a walker and denote by q the probability for the ion to 
occupy the site Q of A neighboring the pore. 

Particles that enter the pore in the low-affinity state 
can exit the system with probability v. With the same 
probability a particle trapped in the pore in the high- 
affinity state can exit A when the pore switches to the 
low-affinity state. We can write the outgoing flux as 
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where we have denoted by r the probability that in the 
stationary state the pore is occupied by a particle when 
it switches from the high to the low-affinity state. 

This one dimensional model can be solved following 
the same idea used in f24]. We find 



N + B + C- y^iN + B + Cy~ 4NC 



and 



2C 



2p 



(1 ~r) N 



(2) 



(3) 



where we have set 

2L^ -3L + 1 



and 



III. 



C = 1 + pAv 



COMPARISON WITH EXPERIMENTAL 
RESULTS 



In 21 it has been demonstrated that the KcsA, a bac- 
terial ion channel of K^^ topology and very well known 
structure (see, for instance, |B] and references therein), is 
a potassium channel. In that paper very precise measure- 
ments of potassium (and some other ion species) currents 
are reported. In this section we try to explain their re- 
sults via our model. 

Since the structure of the KcsA channel is well know, 
in this section will be mainly concerned with the exper- 
imental measures in [2J . At the end a brief analysis of 
the results in [52] and [53], concerning respectively the 
MaxiK and the Shaker channel, will be given. Note that 
the Shaker is a voltage-activated channel of topology. 

A. Potassium current— voltage curves 

In '2r, Figure 2B] the potassium current-voltage curves 
at different concentrations 20, 50, 100, 200, 400, 800 mM 
are shown. We text the reliability of our lattice model 
by exhibiting neat fitting of those experimental data. 

Let / be the measured current, V the voltage applied 
across the cellular membrane and write 



I = Sif 



(1) where / is given in ([!]) with 



V — S\, 



g-a5Vg-Q(l-25)y 



I f,~a(l-2S)V 



(4) 



(5) 



where Si, Sy, S are positive parameters and a is the con- 
stant a = eo/KbT = 0.039 mV"\ with eo = 1.6 x 
10^^ pC the charge of the potassium ion, K^^ the Boltz- 
mann constant, and T = 298.15 °K the temperature. 

As we discussed in the Introduction, the choice of the 
probability v as a function of the voltage V is the only in- 
gredient of the model related to the dynamics of the ions 
inside the filter. The choice ^ is inspired by the model 
in [TMT51 and dates back to the knock-on model in 1261. 



5 





KcsA 


Shaker 


MaxiK 


N 


oUU 


4UU 


zUU 


iUU 


OU 


on 


out) 




/UD 


TO 

1 6 




A nn 
4UU 


1 cn 
ioU 


p;n 


Si (10° pA) 


3.83 


4.91 


6.33 


8.07 


12.0 


16.7 


1.00 


1.03 


1.31 


2.09 


3.03 


7.56 


7.03 


8.88 


P (10-^) 


3.93 


2.07 


1.12 


0.61 


0.32 


0.14 


5.69 


3.91 


2.48 


1.27 


0.46 


0.30 


0.16 


0.08 


s (10"" sec) 


4.17 


3.25 


2.53 


1.98 


1.34 


0.96 


16.0 


15.5 


12.2 


7.66 


5.28 


2.12 


2.77 


1.80 


D (10"^ cmVsec) 


3.62 


4.64 


5.98 


7.62 


11.3 


15.8 


0.88 


0.91 


1.15 


1.83 


2.65 


7.11 


6.60 


8.34 



TABLE I. The parameters Si and p have been measured by fitting the experimental data in [21], [23] (data extracted from 
figure IB therein), and [22] via ^ and ([H]). By fitting the highest concentration set of data we also found Sv = 1.2 x 10^'^, 
S = 0.19, / = 5.5 nm, and L = 1.00 x lO'' for the KcsA, Sv = 5.45 x 10-^ 5 = 0.16, / = 5.3 nm, and L = 1.12 x 10* for the 
Shaker and Sv = 8.91 x 10~^ 5 = 0.16, I = 5.5 nm, and L = 1.00 x 10'' for the MaxiK. The parameters s and D have been 
computed via (|8| and fgf. 



We remark that, as it will be made clear in Section |IIID[ 
the ability of our model to describe the experimental re- 
sults does not depend strictly on this choice. See also the 
comments there in this connection. 

We fit the experimental data for concentration 800 mM 
by using the above formula with N = 800. In this way 
we fix the parameters 

Si, Sv, S, L, and psoo 

where psoo is the low-afHnity probability for concentra- 
tion 800 mM. 

We complete the analysis of the data set according to 
the following scheme. The values of Sj and p for the six 
experimental series are supposed to change, that is to 
say we assume that the probability that the filter is in 
the low-affinity state and the constant Si (we will see in 
Section IIVI that this constant is related to the diffusion 
coefficient in the intracellular region) depend on the ionic 
concentration. The five other experimental series, that is 
to say those measured at concentrations 400, 200, 100, 50, 
and 20 mM, are then fitted by the same equation by 
keeping fixed Sy, S, a, and L and by using p and Si 
as the sole fitting parameters. Results are plotted in 
figure [2j the fitted parameters are reported in table |l] 
The physical meaning of the fitting parameters will be 
discussed in Section HVl 

This fitting scheme is based on a simple remark: we 
assume that the channel behavior, that in our descrip- 
tion is modeled by the function ([5|, does not depend 
on the ion concentration. On the other hand we cannot 
exclude that the ionic diffusion coefficient in the intracel- 
lular region |27| and the typical time spent by the pore 
in the low-affinity state depend on the concentration of 
the ionic species in the cell. 

The soundness of the values that we found for the 



fitting parameters will be discussed in Section jlVj We 
now mention only two facts: to our knowledge there ex- 
ists no experimental measure of the time fraction spent 
by the filter in the low-affinity state. So that we can- 
not say anything about the validity of our estimate for 
p; we can just remark the qualitative agreement with 
PUI Figure 5]. The estimated diffusion coefficient is sev- 
eral orders of magnitude larger than the real (experimen- 
tal) value |17| [5S] . This expected problem is due to the 
one-dimensional character of the model. In Section IIVI 
we shall compute the associated three-dimensional esti- 
mates that will result to be in very good agreement with 
the experimental values. 

Finally we note that with our fitting scheme we also 
are able to reproduce with good accuracy the data in [32] , 
and [23] referring, respectively, to the MaxiK and Shaker 
channel in [23| (see figure |3] and table |l]). Data in figure [s] 
have been reproduced with the authors' permission. The 
data related to the Shaker channel have been extracted 
from figure IB in [23) . 



B. Potassium conductance curves 

In [3T] the permeation behavior of the pore has been 
investigated also by means of the conductance g = I/V. 
We compared the permeation data in [2T] at fixed voltage 
V = 200 mV and concentration varying from 5 mM to 
1600 mM, with the results of our fit of potassium current- 
voltage curves (see Section III A I, as shown in figure |4] by 
an Eadie-Hofstee plot. In the picture the symbols -I- and 
O refer, respectively, to the experimental data and to the 
theoretically computed values, whereas the dotted line is 
an eye guide. 

The matching is good in the whole region where full ex- 
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perimental data sets are known. This graph shows neatly 
the abihty to our model to capture also the dependence 
of the permeation current on the ion intracellular con- 
centration. 

Even if we do not have any physical argument to sup- 
port this choice, it is worth noting that the dotted line 
in the picture, which is just an eye guide, has been in- 
deed obtained by assuming a (slightly sub-linear) power 
low dependence of the low-affinity state probability p on 
the ion concentration and an Hill type behavior for the 
diffusion coefficient. 



C. Other species current— voltage curves 

We test our model against the experimental current- 
voltage curves |21j for ions different from potassium. 
More precisely we consider the curves for NH4" , T1+ , and 
Rb+ at ionic concentration 100 mM. 

To fit these curves we fix the values of S"/, L, and p to 
the ones found previously for potassium at concentration 
100 mM, and we use 5 and Sy as fitting parameters. 
That is to say that we assume that the ionic diffusion 
in the intracellular region and the probability of finding 
the filter in the low-affinity state are independent on the 
ionic species, but depend only on the ionic concentration. 
This is suggested by the measured values of the diffusion 
coefficient, see for instance [251 Table 2] for T1+ and [SDl 
Table 1.1-1] for the other ions. 

With this assumption we are implicitly saying that 
the reduced, with respect to potassium, permeation rates 
typical of the other selected ion species is due to the be- 
havior of such ions inside the selectivity filter. Results are 
plotted in figure [5] the fitted parameters are reported in 
table n 

Deeper investigation on this point should be corrobo- 
rated by additional data. As in the potassium case, data 
sets corresponding to different values of the intracellular 
concentrations are needed. 



D. Modeling the channel 

In the discussion above the choice of the relation between 
the ejection probability and the voltage across the mem- 
brane, see equation ([5]), has been inspired by the physi- 
cal argument proposed in the kinetic model for single-file 
ion channels in p^Hl8] . The function in Q has the same 
form of the channel permeation rate used in that model. 



that is to say we are assuming that the inner part of the 
channel behaves like a single file ion channel. One of the 
key novelties in our work is the coupling between this 
mechanism and the diffusion of the ions in the intracel- 
lular region. 

It is worth noting that the analysis conducted in Sub- 
sections III A |IIIB[ and III C can be repeated using a 
simple power law (having no particular physical meaning) 
for the ejection probability as a function of the external 
voltage. In other words we can assume 

V \T 



(6) 



for modeling the relation between the ejection probability 
and the voltage across the membrane. 

In this case we get similarly good results, that is to 
say we are able to reproduce the experimental data with 
good precision. We do not show the graphs, which are 
similar to those in figure [2j but report the whole set of 
fitted parameters in table |IIl[ 

This remark suggests that, in the framework of our 
model, no particular modeling of the dynamics of ions 
inside the channel is needed in order to get the correct 
behavior of the current -voltage curves. But it is impor- 
tant to remark that different filter models produce differ- 
ent values for the time fraction spent by the selectivity 
filter in the low-affinity state. Indeed, see tables |l] and 
III by using the models ([5| and Q for the dynamics in- 
side the filter, similar values for the diffusion coefficient 
are found, whereas the predictions for the time fraction 
p spent by the selectivity filter in the low-afRnity state 
differ by an order of magnitude. This suggests that the 
current-voltage curves could be a useful instrument to 
predict the time fraction spent by the selectivity filter in 
the low-affinity state once the dynamics of the ions in 
the filter is known or vice-versa. 



IV. DISCUSSION 

In the above section we have seen that the model pro- 
posed in this paper is able to explain in great detail the 
experimental data for the current-voltage curves in the 
KcsA ionic channel. We have shown that the theoret- 
ical predictions fit nicely the measured curves provided 
the parameters defining the model are suitably chosen. 
It is notable that our model explains the dependence of 
the permeation current both on the external voltage (fig- 
ures [2} |3) and [5]) and on the ion concentration in the 
intracellular region (figure |4]). 
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NH4+ 


TI+ 


Rh+ 


5 


0.17 


0.017 


0.33 


Sv 


3.2 X 10"" 


6.8 X 10"* 


9.4 X 10"^ 



TABLE II. Parameters measured by fitting the experimental data for ions different from potassium in KcsA via Q and ([5|. 
Recall also that N = 100, Si = 8.07 x 10^ p = 6.1 x 10"*, s = 1.98 x 10"^* sec, D = 7.62 x 10"^ cm^sec and L = 10*. 





KcsA 


N 


800 


400 


200 


100 


50 


20 


Si (10^ pA) 


4.47 


6.03 


7.71 


9.73 


14.5 


19.6 


V (10-2) 


3.02 


1.77 


1.14 


0.82 


0.45 


0.39 


s (10-** sec) 


3.38 


2.65 


2.07 


1.64 


1.11 


0.82 


D (10-2 cm^/sec) 


4.47 


5.70 


7.29 


9.20 


13.7 


18.5 



TABLE III. The parameters Si and p have been measured by fitting the experimental data in [21] via ([4| and Q. By fitting 
the highest concentration set of data we also found Cv ~ 6.31 x 10*, 7 = 1.19, I = 5.5 nm and L = 1.00 x 10*. The parameters 
s and D have been computed via Q and ([9|. 



In |17| the same experimental data |2lj have been ana- 
lyzed via a kinetic model formerly introduced by the same 
author. In that model the selectivity filter is treated as an 
isolated structure and the intracellular ion concentration 
is an input constant of the model. Two different regimes 
for the filter had to be assumed corresponding, respec- 
tively, to high (400 mM and 800 mM) and low (20 mM, 
50 mM, and 100 mM) potassium concentrations. Two 
different sets of the parameters carachterizing the filter 
behavior (the analogue of Sy and 5) were found. 

In other words the model in [17 predicts that the be- 
havior of the selectivity filter depends on the intracellular 
ion concentration. In our model the selectivity filter is 
coupled with the diffusion of ions in the intracellular re- 
gion; this is indeed the key novelty and the distinguishing 
feature of our approach. This allows to explain the exper- 
imental data referring to all the concentrations with the 
same set of parameters, Sy and 5, for the filter behavior. 

To compare the way in which data are explained by 
the model in |17j and by our model, we first note that 
the electrical dissociation distance 5 (see [151 Fig- 1]) that 
we find is very close to that fitted there. Moreover, we 
note that in [17j in the lower intracellular concentration 
regime, the fitted parameters are such that an higher 
effective concentration (c* in their notation) is seen. In 
our model particle diffusion accounts for this effect. 

We now discuss the reasonableness of the values we 
found for these fitting parameters for the potassium cur- 



tion of the lattice modeling the cytosol. In other words we 
have to associate with the unit length of the lattice model 
a physically reasonable quantity. We imagine to asso- 
ciate with each lattice site a small cubic volume whose 
side length is denoted by I. Recall that the experimental 
concentration c is expressed in niillimolar, that is as num- 
ber of moles per cubic meter (i.e., number of millimoles 
per liter). Recall, also, that N is the number of particles 
in the lattice model and that L is the number of sites in 
A. We then have the following identification 

N 1 

Na¥l 

where Na = 6.022 x 10^^ is the Avogadro number and £ 
has to be expressed in meter. 

Since in our fitting procedure we have chosen N such 
that c/N = 1 mM, we have that 

1 _._io 



rent data for the KcsA ionic channel (see Section III A I 
As a first step we have to give a continuum interpreta- 



=5.50 X 10"*" m = 0.55 nm (7) 

where we have used that L has been fitted to 10^, see the 
caption of table |Tj The above argument provides us with 
a way to associate a length unit measure with our lattice 
model. The result we found for £ is physically reasonable, 
indeed we can assume that each cubic volume associated 
with a site of the lattice can accommodate few (say one 
or two) ions, so that its side length has to be of the order 
of magnitude of the ionic diameter. This value we found 
for £ is then reasonable, since potassium ion, atomic, and 
van der Waals radius are respectively given by 0.13, 0.23, 
and 0.27 nm. 
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As a second step we have to give a continuum interpre- 
tation of the discrete unit of time of the Markov Chain. 
In other words we have to associate with the unit of time 
a physically reasonable quantity s. We consider the evo- 
lution of our model up to the time n and note that the 
true current measured in the experiments can be identi- 
fied as 

^ eg X (number of particles which exited A) cq ^ 
s X (number of iterations) s 

where bq — 1.6 x 10^^ pC is the charge of a potassium 
ion and / is the theoretical model outgoing flux. We 
then identify cq/s with the parameter Sj introduced in 
equation Q; that is to say we write 



eo 
Si 



(8) 



and we can use the estimated Sj listed in the flrst row 
of table |l] to compute the values of s listed in the same 
table. 

We have no direct clue to establish if these values for 
s are reasonable or not. But starting from this values we 
can estimate the diffusion coefficient of potassium ions 
in the cytosol. Indeed, see Appendix [B] the intracellular 
ion diffusion coefficient is related to the other parameters 



by the equation (Bl), which yields 

£2 



D = 



2s 



(9) 



By using this equation we compute the diffusion coeffi- 
cient for the potassium. Results have been reported in 
the last row of table H] 

The order of magnitude we found for the diffusion co- 
efficient of potassium ions is 10"^ cm^/sec. This result 
is several (say three) order of magnitude larger than the 
real (experimental) value |28) . This problem was indeed 
quite expected, since we modeled the diffusion of ions in 
cytosol with a one-dimensional lattice. As we have al- 
ready remarked, the one-dimensional choice is motivated 
by the possibility to write explicitly the solution of the 
probabilistic model. Explicit formulas are quite neces- 
sary to perform the extended analysis discussed in Sec- 



tion III As a future work we are now planning a Monte 
Carlo assisted analysis for the three-dimensional version 
of our model. 

In any case we prove, now, that the estimates we found 
in this one-dimensional case for the diffusion coefficient 
are quite reasonable. Consider a simple symmetric ran- 
dom walk with unitary time on a three-dimensional cu- 
bic lattice with L sites and spacing m. We consider the 



side length L^^^ in order to ensure that the "real" vol- 
ume associated with the three dimensional model is equal 
to that associated with the one dimensional one. By a 
classical argument similar to the one developed in the 
Appendix |B] it is proven that the squared mean distance 
walked up to time t is QD^t where is the diffusion 
coefficient. This implies that the typical time needed by 
the walker, started at random in the lattice, to reach a 
particular point of the boundary (say the pore) is of order 

^ X X 6 i^/^ 2 

64 6D3 ^ ' 

Indeed, the first term is the square of the average dis- 
tance of a point inside a cube of side length mL^^^ from 
the boundary of the cube itself and the third one is the 
number of times the walker has to visit the boundary 
before touching the pore (number of points on the lat- 
eral surface of the cube). Note that the product between 
the first and the second term is an estimate for the time 
needed to touch the boundary of the cube for a walker 
started at random in the cube itself. 

In our former analysis, in computing fluxes, we in- 
deed evaluated this time in the framework of our one- 
dimensional model. By repeating the same argument we 
can say that this time was estimated as 



X X 2 
2Di 



By equating the two expressions found above we find 
1 



-DiL 



-2/3 



0.00054 X L>i 



where we used L = lO** (see the caption of table |l|. 
By the equation above it follows that to the diffusion 
constant estimated via the one-dimensional model rang- 
ing in 3.62 - 15.8 x 10^^ cm-^/sec (see table |l| it cor- 
responds a "true" diffusion coefficient ranging in 1.95 - 
8.53 X 10~^ cm^/sec, which is very close to the experi- 
mental value for the potassium diffusion coefficient [28] . 



V. CONCLUSIONS 

In this paper we have studied the effect of intracellu- 
lar ion diffusion on ionic currents permeating through 
the cell membrane ionic channels. Ion channels share 
the following common properties: the ion flux is rapid, 
the channel is selective, and its functions are regulated 
by a gating mechanism. Although a formidable effort, 
with absolutely remarkable results, has been performed 
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recently, a complete understanding of ionic channel be- 
havior is still lacking. 

Thanks to the patch clamp technique very precise mea- 
surements of the true current across a single channels are 
available. In this paper we have proposed a model which 
is able to provide a full explanation of the KscA current- 
voltage experimental curves j21i by taking into account 
the following three effects: ion diffusion in the intracel- 
lular region, dynamics of the filter, and dynamics of the 
ions inside the filter. 

In particular our model points out the role played by 
intracellular diffusion on ion permeation in compensat- 
ing the ion depletion in the region close to the pore. 
The model is able to explain the dependence of the ex- 
perimetal data referring to the true ionic currents both on 
the external voltage and on the intracellular potassium 
concentration. 

Is is also notable that all the fitting parameters have 
a clear qualitative and quantitative physical interpreta- 
tion. Moreover, we noted that the estimate we got for 
the typical time spent by the filter in the low-affinity 
state is strictly related to the dynamics of ions inside the 
filter. This suggests that our approach, corroborated by 
new experimental data, should be able to shed some light 
also on the dynamics of the ions inside the filter. 

Appendix A: Detailed definition of tiie model 

We consider an integer time variable n. We set n = 
and choose at random with uniform probability 1/L the 
position of the N particles. We then repeat the following 
steps: 

1 . set n equal n + 1; 

2. select at random the state of the pore: choose the 
low-affinity state with probability p and the high- 
affinity one with probability 1 — p; 

(a) if the pore is in the low-affinity state and it is 
occupied by a particle, the particle is released 
with the following rule: it jumps with prob- 
ability 1 — V to the site of A neighboring the 
pore or, with probability v, it exits the system; 

(b) if a particle exits the system, a particle of the 
same species is put at random with uniform 
probability 1/i on one of the L sites in A; 

3. the position of each particle on the lattice is up- 
dated following the rules defined in Section |llj 

(a) if a particle enters the pore and the pore is in 
the low-affinity state, the particle is immedi- 
ately released by the pore with the following 



rule: it jumps with probability 1 — v to the site 
of A neighboring the pore or, with probability 
V, it exits the system; 

(b) if a particle exits the system, a particle of the 
same species is put at random with uniform 
probability 1/i on one of the L sites in A. 

Appendix B: Diffusion coefficient 

Consider a one-dimensional walker on iZ and denote its 
position at time t G sZ+ by 

Ste{...,-£,0,+i,...} 

The random variable St is the sum of t/s (note that t/s 
is a positive integer) independent identically distributed 
random variables Yi E {—£, +^} such that 

p(y, = +i) = p(r, = -i) = l 

A straightforward computation yields 

E[Y,] = and E[{Y, - E[Y,])^] = 

for the average and the variance of the random variable 
Yi respectively. By independence we also get 

E[St/s] - and £[(5*/, - E[St,,]f] - U'' 

By the Central Limit Theorem (see, for instance, [5T1 
Theorem (4) in Section 5.10]) we have that for t large 
the probability density of the random variable St/s is 
very well approximated by the Gaussian 

, ^ exp{-2;V[2£2(^/^)]| 

for a; €E M. By comparing this result with the continuum 
description given by a one-dimensional diffusion equa- 
tion u — Du" with diffusion coefficient D we get the 
identification 

2Ds = (Bl) 
allowing to compare lattice and continuum results. 
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High-affinity: Low-affinity: 

o— e— ^ o^-e — • o >• 

(free) (occupied) 



FIG. 1. Rule for the site close to the pore. The black solid circle denotes the pore, the open circles denote sites of the lattice, 
the black arrows denotes jumps that arc performed with probability 1/2, the black lines denote jumps that cannot be performed 
(probability zero). On the left the behavior of an ion close to the pore in the high-affinity state is depicted: if the pore is 
free, the particle jumps with uniform probability 1/2 to one of the 2 nearest neighboring sites; if the pore is occupied, the 
particle jumps with uniform probability 1/2 to the nearest neighboring site in the lattice (it cannot enter the pore) and with 
probability 1/2 the particle does not leave the site. On the right the behavior of an ion close to the pore in the low-aiiinity 
state is depicted: the particle jumps with uniform probability 1 /2 to one of the 2 nearest neighboring sites; If the pore is in the 
low-afHnity state, the ions behave with the same rule as that for an ion faced to the free pore in the high-affinity state. It is 
worth remarking that the pore in the low-affinity state cannot be occupied. 
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FIG. 2. Comparison between experimental KcsA potassium current-voltage curves in |5T] (symbols) and theoretical prediction 
(solid lines). The symbols Q, M, □, *, x, and + refer respectively to the potassium concentrations 20, 50, 100, 200, 400, 800 mM. 
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FIG. 3. Comparison between experimental MaxiK potassium current -voltage curves in 22 (symbols) and theoretical prediction 
(solid lines). The symbols *, x, and + refer respectively to the potassium concentrations 50, 150,400 mM. 
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FIG. 4. The Eadie-Hofstee plot for the potassium permeation current at external voltage 200 mV. The symbols + and O 
refere, respectively, to the experimental data and to the theoretically computed values. The dotted line is an eye guide. 
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FIG. 5. Comparison between experimental current -voltage curves [21] (symbols) and theoretical prediction (solid lines). The 
symbols *, +, and x refer respectively to Rb"*", NHj, and T1+ at concentration 100 mM. 



